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Abstract 

This paper considers passing from the usual R d model of absolute space to at the level of relational particle 
models. Both approaches’ d = 1 cases are rather simpler than their d > 2 cases, with N particles in S 1 admitting 
a straightforward reduction with shape space T A_1 . The § 2 and § 3 cases - observed skies and the simplest closed 
GR cosmologies respectively - are also considered, the latter in the contexts of both static and dynamical radius of 
the model universe. The space of relational triangles on § 2 is hyperbolic 3-space H 3 . Overall, by passing to a closed 
underlying absolute space, and then to dynamical notion of space, we close some of the modelling gaps between 
relational particle models and geometrodynamics or its inhomogeneous perturbative regime of interest in cosmology. 
Quantum counterparts are also outlined, for use as model arenas of quantum cosmology. These models are useful in 
further considerations of both classical and quantum background independence. 


1 Introduction 

Newtonian Mechanics, and the Newtonian Paradigm for Physics more generally, are based on absolute space and time. 
The immovable external character of these led to these being opposed by relationalists, in particular Leibniz [1] and 
Mach [2]. On the other hand, the Newtonian paradigm of Physics sufficed to explain the observations of nature until 
the end of the 19th century. Furthermore, a satisfactory relational alternative to the foundations of Mechanics was not 
found until 1982 by Barbour and Bertotti [3]. Of course, this long post-dates 1) the advent of GR and 2) its dynamical 
reformulation [4], by which incorporated a number of relational features into Physics: closed universes, dynamical rather 
than background geometry. Yet, it turns out that relational particle mechanics (RPMs) [3, 5, 6, 7, 8] are useful model 
arenas as regards understanding many a feature of General Relativity (GR) as a dynamical system [4, 7], of GR’s 
configuration spaces [9, 10, 8]. And of analyzing which aspects of Background Independence GR possesses, difficulties 
with which then become facets of the notorious Problem of Time in Quantum Gravity [11, 12, 13]. Indeed, GR admits 
a geometrodynamical formulation [4] that can be cast in Machian form [12]. This in the sense of obeying postulates of 
Temporal and Configurational Relationalism (Sec 2) which can already be set up for RPMs. Thus the study of RPMs 
turns out to have further foundational value (also the above two postulates are also a subset of those of Background 
Independence [11]). 

Each of [8] and the current Article extend the range of such relational models in different ways. This is to better cover 
a number of the significant gaps between RPMs and GR. The present paper addresses this by use of curved geometry, 
of closed-universe models, and of dynamical geometry: a range of further models lying in between the two. Firstly (Sec 
3), I consider relational theory of point particles where S d plays the role of absolute space in place of the habitual 
of Mechanics. These could be the very simple S 1 [7], S 2 either as a model or as the observed sky [14, 15, 7], or the S 3 
that is one of the simplest closed models of space in GR cosmology. In Sec 4, I consider the geometry of the spaces 
of points on a circle and of spherical triangles, alongside Shape Statistics thereupon [14, 15]. Shape Statistics of point 
configurations on the sky is of value to perturbatively inhomogeneous Cosmology [16, 25, 10], and means of probing 
further tests, concepts with interpreting the CMB and the distribution of galaxies. Secondly, I consider relational theory 
on expanding spheres (Sec 5). I end by outlining quantum version of spherical RPM in the sense of Dirac quantization 
in Sec 6. 


2 Outline of the Relational Approach to Physics 

The first structural element in the approach is configuration space C| [18]: the space of generalized configurations Q A 
for a physical system. This paper considers the case of point configurations. Here, an incipient C| is 

c\(N,d)=xf =i a(d) (i) 

for H(d) the corresponding notion of absolute space. 

Configurational Relationalism is then that 

i) no extraneous configurational structures - spatial or internal-spatial metric geometry variables that are fixed-background 
rather than dynamical - are to be included in the formulation. 

ii) Physical formulations in general involve not only a C| but also a 0 of transformations acting upon CJ that are taken 
to be physically redundant. 
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For point particle models, ii) alone is realized [i) is also realized for GR, where the configuration space metric is a 
function of the 3-metric configuration]. In some cases, 0 can be implemented directly, due to enough 0-invariant objects 
being available. More widely, 0 can be implemented indirectly, whether by A) Best Matching [3], which involves 0 
corrections to the theory’s changes of redundant configuration, dqj (see Sec 3 for examples). Or B) by the more general 
‘0-act 0-all’ method (see e.g. [8]). 

On the other hand, Temporal Relationalism - the Leibnizian stance that that there is no time at the primary level for 
the universe as a whole - can be formulated as follows. 

i) The action is not to include any extraneous times - such as Newtonian time - or extraneous time-like variables, such 
as the lapse function [4] in the geometrodynamical formulation of GR. 

ii) Time is not to be smuggled into the action in the guise of a label either. 

Here, ii) is most optimally implemented by a geometrical action, which is dual to parametrization irrelevant action, which 
upon parametrization would become a reparametrization invariant action [7, 19]. The particular cases of geometrical 
actions in the present Article are all Jacobi-type action [18], 

S = V2 J d sVE - V , (2) 

corresponding to a Riemannian notion of geometry. A primary constraint must then follow from these actions due 
to an insight of Dirac’s [20]. In the case of Mechanics, this gives an equation which is usually interpreted as an 
energy conservation equation. However in the present context this is to be reinterpreted as an equation of time [21]. 
Indeed rearranging it gives an expression for emergent Machian time: a concrete realization resolution of primary-level 
timelessness by Mach’s ‘time is to be abstracted from change’. In the case of metric Scale-and-Shape Mechanics, this 
emergent time amounts to a relational recovery of a quantity which is more usually regarded as Newtonian time. Note 
that more generally handling Mechanics models in a temporally relational manner requires a modified version of the 
Principles of Dynamics as laid out in [19]. 

Another useful notion is relational nontriviality, which requires at least 2 degrees of freedom, so that the value of one 
can be expressed in terms of the other (rather than in terms of a no longer a priori existent external time parameter). 
This leads to the notion of a smallest relationally nontrivial unit for each relational theory [8]. This is concurrently 
the smallest relational universe model, the smallest relationally nontrivial subsystem, and the smallest Shape Statistics 
sampling probe. An archetype for this is the relational triangle in 2- and 3 -d metric scale-and-shape and pure-shape 
RPMs, as features in both Barbour’s familiar demonstration of Best Matching with two wooden triangles and in Kendall’s 
sampling by triples of points in 2 -d. Furthermore, the geometry of the shape space of these probes is geometrically well- 
understood: it is S 2 . (Or a piece thereof: Kendall’s spherical blackboard [22, 14]. See [10] for discussion of various pieces 
of this that are relevant in different modelling contexts.) This S 2 geometry then enters many subsequent calculations 
for the corresponding Shape Mechanics and Shape Statistics [15, 7, 25]. 

3 RPMs for particles on spheres 

3.1 General case 

Here absolute space is modelled by Cl = S p , 

q(N,d)=x* =1 sr , (3) 

and a natural and nontrivial choice of 0 is 

Isom{ S p ) = Rot{p + 1) = SO{p + 1) . (4) 

The relational action for this is built out of a potential function V = V{p~<) - i.e. a function of relative p-sphere angles 
alone - and a curvilinear rendition of Best Matching based upon SO(p + 1) auxiliaries S^: (2) with 

d « 2 = 5Zf=i dss2 ' ( 5 ) 

Then the quadratic constraint 

£ : = E 2 + V ^~<) = E (6) 

arises as a primary constraint, for 

N pIqJ := TT P- 1 , sin ~ 2 e mI S pq 5 IJ (7) 

X X m = 1 
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the inverse of the standard p-sphere metric, 

MpIqJ := ^T^OmlSpqSlJ . ( 8 ) 

(i 9 mI here are N copies of the p-sphere angle coordinates and E is total energy of the model universe.) 

Also the linear constraint M t j = 0 arises as a secondary constraint from variation with respect to the Shy. 

Moreover (6) can indeed be rearranged to give an expression for classical Machian emergent time: 

^em(JBB) = E SeSO(p+1) ds / ^2{E-V(^<)}^ . (9) 

ds here depends implicitly on the SO(p + 1) auxiliaries via (5), whereas E denotes ‘extremization over’. 

The minimal relational nontrivial units have 2 < dim(C|) — dim(fi) = Np — {p + l}p/2 = p{2IV — p — l}/2, giving 3 
points on the circle for p = 1, the spherical triangle for p = 2, and the 3-spherical triangle for p = 3. 


3.2 Specifics of S 1 ‘circle lines’ 

In this case, absolute space is modelled by Cl = S 1 , 

q = xf =i s 1 , (to) 

and a natural and nontrivial choice of is 

/som(§ 1 ) = Rot{ 2) = 50(2) = 17(1) . (11) 

Extending the nomenclature of [23], who considered particles on a line as ‘metro line configurations’, the current theories’ 
configurations are topologically circular lines, such as indeed the London Underground’s circle line. The spaces of these 
configurations then form IV-stop metroland for the ffi. case and N-stop circle-line-land for the new S 1 case. 


The relational action for IV-stop circle-line-land is then built out of V = ^(relative angles (j) 1 
Matching 

ddi := dd 1 - da : 

(2) with 

ds2 = • 


The conjugate momenta are then 


These obey 


p , = ME-VW-f )} , 

as 

£:= Ef=i pV/z+v ^ 1 -t J ) = E 


and the Best 

( 12 ) 

(13) 

(14) 

(15) 


as a primary constraint, and c = 0 as a secondary constraint from variation with respect to the single Rot{2) auxiliary, 
a. 


For these models, it is immediate to solve c = 0 as an equation for a, giving 

da = ” = 1 d( t ,I / N : = W) : ( 16 ) 

an averaged object. Then substituting back into the kinetic arc element, 

d?2 = EL d ^ 2 ( 1? ) 

for contrast objects defined according to O := O — (O). There are n = N — 1 independent such, amounting to losing one 
S 1 in the quotienting. In some ways, this is the circular analogue of passing to centre of mass frame in quotienting out 
translations. 

Finally note that (15) can be rearranged to give 

< em = E aeflot(2) (y ds/yjtiE-VW-V)}} = j ds/^2 {E-VW-p)} (18) 

for ds depending implicitly on a through (23). The last more explicit line is obtained through the above reduction. 
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Figure 1: A circle line configuration, and a minimal relationally nontrivial spherical triangle configuration. 


3.3 Specifics of § 2 ‘constellations’ 


Now, absolute space is modelled by d = S 2 , 


q = xf =i s 2 , 

(19) 

and a natural and nontrivial choice of is 


Isom( S 2 ) = Rot( 3) = SO(3) . 

(20) 


It is particularly appropriate then to term the point particle configurations in this case constellations, the corresponding 
space of which is constellationland. 


The corresponding rendition of Best Matching corrections now take a curvilinear form: 

d bO 1 := d 0 1 + cos (j ) 1 AB 2 — sin </> 7 AB\ , 

Agf /) 1 := Ac /) 1 + AB 3 — cot0 7 {sin</> 7 AB 2 + cos (j ) 1 ABi} . 
The relational action is (2) with V = V (<) alone and 

ds 2 = ^2 ^ =1 {d B 9 2 + sin 2 6*/dij(/> 2 } . 

Then the conjugate momenta are 


«_ VHe-v«)} a aI 


Pi = 


, al <t> ■ 2 Q y/2{E-V(<)} , Al 

A B 9 , p} = sm 0/ -- - - A B (j> . 

As ~ As ~ 


These obey 


sin ~ 20/ Pi*}/ 2 + V (<) = E 

as a primary constraint, and 

c = ^2 —sin — cos ^cot dipf , cos — sin c /) 1 cot dipf , pf 

as a secondary constraint from variation with respect to B. 

Furthermore, (25) can be rearranged to give the emergent Machian time 

< em = Eb 6S o(3) As/^/2{E-V{<)} 
for ds depending implicitly on as per (23). 


= 0 


( 21 ) 

( 22 ) 

(23) 

(24) 

(25) 

(26) 


(27) 


3.4 2-d with other metrics 

As a useful digression at this point, the general ellipsoidal line element is [e.g. using x = (a sin 9 sin <fi , b sin0 cos <j>, c cos 9) 
in a fiducial embedding flat ]R 3 ] 

ds 2 iiip Soi d = {cos 2 0{a 2 sin 2 </> + 6 2 cos 2 ^} + c 2 sin 2 0}d0 2 

+2{a 2 — 6 2 }sin 9 cos 9 sin</>cos (f>A9A<p + sin 2 0{a 2 cos 2 0 + b 2 sm 2 9}A(f> 2 . (28) 

Now for the spheroidal special case a = b, the line element collapses to 

ds 2 pheroid = {“ 2 cos 2 6»7 + c 2 sin 2 0/}d0 2 + a 2 sin 2 0/d((> 2 , (29) 
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so (j) is a cyclic coordinate, corresponding to an axial symmetry SO( 2). In this case (and for a = b so as to not fall 
back to the previous Sec’s very special case), 

/som((S 2 , dSg pheroid )) = 50(2) (30) 

is available as a 0 of physically meaningless transformations. Then take C| as per (19) with each S 2 equipped with this 
metric. This amounts to summing over N labelled copies of this metric, alongside applying the Best Matching correction 
d b^ 1 '■= d<7 + d B 3 . This makes sense in conjunction with the 50(2) respecting form V = P(# 7 alone). 

On the other hand, for the generic ellipsoidal case - a, b,c all distinct - 

/som((S 2 ,d4iip SO id» = ^ , (31) 

so no continuous 0 is available in this case. This is indicatory of the generic situation for RPMs on fixed manifolds: most 
fixed manifolds have enough structure that there is no scope for any Configurational Relationalism at all. Moreover, 
this is a limitation of RPM modelling which the GR situation itself does not succumb to. 

The above examples indicate that there can be multiple metrics on a manifold, which exhibit, moreover, a diversity of 
isometry groups. 

3.5 3-d closed universe particle model 

Here absolute space is modelled by H = S 3 , 

q=Xf =i S 3 , (32) 

and a natural and nontrivial choice of 0 is 

Isom( S 3 ) = 50(4) . (33) 

The relational action is then built out of V = P(3-<) and another use of Best Matching in curvilinear form. The action 
here is (2) with 

ds 2 = ^ ^_ 1 {d sdj + sin 2 0 7 ds0 2 + sin 2 0 7 sin 2 (/> 7 d,s<4} (34) 

now involving two 3-vectors’ worth of auxiliaries .SV^-i: 5 and B 1 , due to the underlying ‘accidental relation’ 

50(4) 3* 50(3) x 50(3) . (35) 


The ensuing constraints are 

£ : = ^2f =1 iP e i 2 +sin- 2 6 ,/ pf 2 sin- 2 (9 / sin- 2 4l5/ 2 }/2 + y(3-<) = E (36) 

as a primary constraint, and the 50(4) constraint Mij = 0 reducing to the schematic form c = 0 and d = 0, as a 
secondary constraint from variation with respect to 13 and id. 

Also now 

t cm = SO (3) X SO (3) (yj &s/\j2{E - P(<)}^ (37) 

for ds depending implicitly on 5 and id through (34). 

Finally note that subgroups of SO(d+ 1) present further options, and that affine and conformal extensions of Relation¬ 
alism - considered in [8] for - are also available for S d . 

4 Shape space geometry and Shape Statistics 

The general shape space for N particles on S d is 

C(N, d) := X ^ =i S d /SO{d + 1) (38) 

In the 1 -d case, 

circle-line-land C(N, 1) := X^3750(2) = X^S 1 = T n (39) 

and 

Isom(C(N, 1)) = Isom{ TP) = X” =i 17(1) . (40) 
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Figure 2: a) Clumping on a line, b) Set-up for angular Shape Statistics in R 2 space, c) Clumping on a circle, d) Set-up for angular Shape 
Statistics in § 2 space; great sphere geodesics have taken the place of lines, and ‘decorated’ refers to non-principal strata. 


On the other hand, in the 2 -d case, constellationland C(N,2) := S 2 /SO( 3). The smallest relationally nontrivial 

such is spherical triangleland C(3, 2) := X 3 S 4 /50(3), whose principal stratum [10] is of dimension 2 x 3 — 3 = 3. 

Furthermore, this is geometrically [15] the hyperbolic space 

H 3 := {x € M 4 |(||:r|| 2 = -l,x 0 > 0} , (41) 

(here M 4 is mathematically the same as 4 -d Minkowski spacetime), with standard line element 


d» = -*i!L 

1 — \\x r 


Then 


Isom{ H 3 ) = 2 -d Mobius group 
the proper orthochronous Lorentz group so familiar from SR. 


50+(3,1) 


(42) 

(43) 


In contrast, the metric shape space 5(N, d) on R d with 5(N, 1) = S n_1 and S(N, 2) = CP” -1 [15, 7] . This shows some 
capacity for arbitrary -N shape spaces being geometrically standard (though this does not carry over to 3-d). Some 
remaining questions then concern the topological and geometrical identity of C(N, 2) for N > 3, and of the probe space 
C(N, 3). Reduction of the Mechanics action can on occasion point to these features [6, 7[. Though metric shapes on S 2 
parallel metric shapes on R 3 [both involve 50(3)] in only having locally attainable reduction, and metric shapes on § 3 
are also expected to follow suit in this regard. 


As regards clumping statistics on circles, Roach’s [24] black and white square procession approach on R 1 can be extended 
to this case. This extension involves truncating the geometric progression of black-square chain-lengths and averaging 
over choice of starting-point for the procession on the circle is performed. 


As an alternative, one could consider Shape Statistics based on Geometrical Probability on S(N,1) = § n : N -stop 
metroland [23], with circle-line-land counterpart of Shape Statistics based on Geometrical Probability on C(N, 1) = T": 
IV-stop circle-line-land. 


Shape Statistics on circles additionally involves relative angle information (the circle clumping case’s angle is but an 
extrinsic notion). This parallels how R 1 has just clumping (ratio) information whereas R 2 has both ratio and angle 
information, thus necessitating angle-based Shape Statistics like Kendall’s [22, 15]. The § 2 angular Shape Statistics 
can be use e.g. to investigate quasar alignment in the observed sky [26]. More generally, shape space geometry and 
Shape Statistics provide the means of devising further and mathematically better-founded pattern tests along the lines 
of ‘circles in the sky’ or ‘multiple image identification’ [27]. 


5 Minisuperspace-RPM coupled models 

One could furthermore entertain the notion of space £ itself being dynamical: a GR-type feature. E.g. £ could undergo 
whichever of expansion, anisotropic change, and inhomogeneous change. A further source of variety in these cases is 
whether to model particles on a dynamical £, fields on a dynamical £ or just the dynamics of £ itself. Consider for 
now an isotropic GR model coupled to spherical RPM. These are in essence a combined minisuperspace-RPM model, 
with 

S = 2 J \j ^ — a 2 —^^\/2Ka{Aa 2 — 3) + V(3-<, a) (44) 
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Here n is the gravitational part of the action’s dimensionful prefactor, a is the cosmological scalefactor and A is the 
cosmological constant. 


The conjugate momenta are then 


yj2na{A.a 2 — 3) + V(3-<, a) na , a \f2na\Aa 2 — 3) + V (3-<, a) 2 - T 

- Tg -yck, Pli = - - - aMujjdP . 


These obey the primary constraint 


< I r 
+ — 
na 


IxfUSS 2 


2a 2 


= 2na{3 — Aa 2 } — V (3-<, a) . 


Finally, 


t em = 


—nd a 2 + 2a 2 ds| 


s 3 


, iem 
c 0 — 


l 


a ds§3 1 


V(3-<) 


4{2«a 2 {3 — Aa 2 } — H(3-<, a) u 2\[2 J Aa 2 — 3 da 128 yfn J a 4 {3 — Aa 2 } 2 

(compare the expansions of classical Machian emergent time in [28]). 


(45) 


(46) 


(47) 


Note 1) As regards the particles breaking the homogeneity, in any case we want the inhomogeneity to be small for the 
cosmological applications in hand, so the spherical RPM’s are to be treated as perturbations about minisuperspace. 
Note 2) A limitation in this analogy is that GR’s inhomogeneous perturbations about minisuperspace have a further 
scalar-vector-tensor (SVT) split of modes. 

Note 3) On the other hand, a virtue of this Sec’s model arena is the availability of further pattern analysis for it, of use 
in analysis of cosmological structure formation. This remains a distant prospect for actual inhomogenous GR. 


6 Quantum spherical RPMs 

Dirac quantization is straightforward. Here the standard spherical kinematical quantization [29] will do. 


In the 1 -d case, this gives 


Cd! = 


h ST N 

i 2.^1 1 = 1 dtf) 1 


= 0 


- meaning that 4/ = 4/(((/ — cf> J ) alone - and the main wave equation 


(48) 


2d! = 


h \ ' N d 2 # 
2 / j J_ \ dQ 12 


+ V(<1> 1 


cj> J )d! = Ed! . 


(49) 


On the other hand, the 2-d case gives 


2 * GEL 


—sin 


iimL 
0 de 1 


COS ((/cot 0/ 


cos 


lIEE 
0 de 1 


sin ((/cot 9 1 


d 

dcfi 1 



(50) 


meaning that \F = 4/(<) alone, alongside the main wave equation 


^ = -§£E^+n<)* = £*- ( 51 ) 

The 3-d case then necessitates 6 components’ worth of 

Mijd> = 0 : 241 = 0 and c/d! = 0 , (52) 

meaning that \F = 4/(3—<) alone, and the main wave equation 

2d! = -f + = Ed! . (53) 

Finally the 3-d case coupled to a GR-type dynamical scalefactor gives (52) alongside the new main wave equation 



^2 E/Li^s 3 2a 2 j 1 H na ^ + 2Ka l A “ 2 " 3 1^ + ^(3-<,a)4 r - 0 • ( 54 ) 

(This makes use of the conformal operator ordering [30].) This would be solved to zeroth order in terms of a and then 
subjected to inhomogeneous RPM particle perturbations. 
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